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The  Heisenberg  exchange  Hamiltonian  has  been  solved  exactly  for  a  cubic  array 
of  eight  spins  (each  with  spin  1/2).  Both  energy  eigenvalues  and  thermodynamic 
functions  have  been  calculated.  A  Curie  and  a  Neel  temperature  can  be  defined 
and  their  values  determined  as  a  function  of  the  strength  of  first,  second  and 
third  neighbor  interactions.  For  some  values  of  the  exchange  constants,  “spiral” 
antiferromagnetic  states  exist. 
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FERRO-  AND  ANTIFERROMAGNETISM 
IN  A  CUBIC  CLUSTER  OF  SPINS 


I.  INTRODUCTION 

The  exact  three-dimensional  solution  of  the  Heisenberg  model  for  ferromagnetism  is  not 

tractable  for  an  infinite  lattice.  There  are  various  approximate  solutions,  some  of  which  in- 
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volve  the  exact  solution  for  some  small  cluster  of  spins.  '  In  addition,  a  number  of  solutions 

3 

exist  in  the  literature  for  a  small  number  of  spins  in  various  configurations.  This  report  gives 
the  exact  solution  to  the  Heisenberg  Hamiltonian  for  eight  spins,  each  with  spin  i/z,  located  on 
the  corners  of  a  cube.  Both  the  eigenvalues  and  the  thermodynamic  functions  have  been  calculated. 

This  simple  cluster  shows  many  of  the  features  of  an  infinite  magnetic  material.  In  partic¬ 
ular,  a  ferromagnetic  Curie  temperature  T^  and  an  antiferromagnetic  Neel  temperature  Tj^  can 
be  defined,  and  their  values  computed  as  a  function  of  the  strength  of  second  and  third  neighbor 
interactions.  It  is  shown  that,  for  some  values  of  the  interactions,  "spiral"  antiferromagnetic 
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states  exist.  It  is  also  shown  that  there  are  regions  in  which,  as  the  temperature  is  lowered, 
ferromagnetic  ordering  begins,  but  at  still  lower  temperatures,  the  system  drops  into  the  anti¬ 
ferromagnetic  singlet  state. 

The  approximate  solution  to  the  infinite  simple  cubic  lattice,  using  the  results  contained 
herein,  will  be  published  in  a  subsequent  report. 

n.  CALCULATION 

The  Heisenberg  Hamiltonian  is  written 

5C  =  —  2J[€  +  xp  +  yv]  (1) 

in  which  J  is  the  exchange  constant,  x  and  y  are  proportional  to  the  strength  of  next  nearest 
and  third  nearest  neighbor  interactions,  and 

e  =S,  •  S^  +  S^-  S3+S3  .  S4  +  S4  •  +S5  -  +S^  •  S^+S^.  Sg  +  Sg  .S5 

+  ■  Sg  +  •  Sg  +  S3  ■  S.^  +  •  Sg 

=  ^1  •  ^3  +  ^2  •  +  S5  •  s.^  +  Sg  •  §g  +  •  Sg  +  §2  •  S.^  +  S3  •  Sg  +  S4  •  Sg 

+  •  ?g  +  ■  ^5  +  ^3  ■  ^6  ^4  ■  ^7 

V  =  S^  •  ?.^  +  §2  •  Sg  +  ?3  •  ?5  + 
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Fig.  2.  Energy  spectrum  assuming  only  nearest  neighbor  interactions, 
X  =  y  =  0.  The  spectrum  is  arranged  according  to  the  value  of  the  total 
spin.  The  N^el  state  energy  is  indicated  even  though  it  is  not  a  proper 
eigenstate  of  the  system. 
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The  spin  operator  for  the  site  follows  the  numbering  of  the  sites  shown  in  Fig.  1, 

The  eigenfunctions  of  Eq.  (1)  can  be  classified  according  to  the  total  spin  S,  the  z-component 
of  the  total  spin  Mg,  and  the  irreducible  representation  of  the  cubic  group.  Table  I  lists  the  sym¬ 
metry  types  and  energy  eigenvalues  for  the  cubic  cluster.^  Appendix  A  tabulates  the  wave  func¬ 
tions.  Appendix  B  gives  the  matrix  elements  of  e  +  xp  +  yi'  for  the  wave  functions  of  Appendix  A. 
The  notation  for  the  representations  is  such  that  a  vector  transforms  as  T^.  The  spectrum  for 
X  =  y  =  0  is  shown  in  Fig.  2. 

For  a  given  value  of  the  total  spin  S,  the  center  of  gravity  of  the  energy  levels  is  given 
by  the  sum  rule 


Ecg=-2 


Z  J-- 

^  11 

Li<i  . 


[S(S  +  1)  -  (3/4)  N]/[N(N-1)] 


(3) 


in  which  N  is  the  number  of  spins  (each  with  spin  l/2)  and  J.^  is  the  exchange  constant  between 
sites  i  and  j.  Equation  (3)  represents  a  convenient  check  of  the  calculated  energy  eigenvalues. 
The  Heisenberg  approximation  for  ferromagnetism^  corresponds  to  calculating  the  partition  func¬ 
tion,  using  Eq.  (3)  for  the  energy  eigenvalues. 

The  partition  function  for  the  cubic  cluster  is  calculated  by  using  the  eigenfunctions  listed 
in  Table  I.  From  the  partition  function,  the  susceptibility  and  heat  capacity  for  the  system  are 
calculated.  The  susceptibility  x  given  by  the  Curie  law  is  written  as 


X  =  NS(S  +  1)  (4) 

in  which  N  is  the  number  of  spins,  S  is  the  spin  quantum  number  g  =  2  and  |Xg  is  the  Bohr  mag¬ 
neton.  The  cubic  cluster  obeys  the  Curie  law  with  suitable  identification  of  N  and  S.  For  high 
temperatures,  2j/kT  «  1,  the  necessary  interpretation  of  N  and  S  is 

NS(S  +  1)  =  6  , 

and  at  low  temperatures,  2j/kT  »  1,  the  interpretation  is  changed  to  N  =  1  with 
NS(S  +  1)  =  20  J  >  0 


and 


NS(S  +  1)  =  0  J  <  0  . 

A  plot  of  the  susceptibility  vs  l/T  for  x  =  y  =  0  is  shown  in  Fig.  3  for  both  the  ferromagnetic  and 
antiferromagnetic  states.  For  convenience,  the  plot  is  presented  in  dimensionless  form  in  terms 
of  «x  vs  2j/kT,  with  a  =  bj/g^p^.  For  the  ferromagnetic  case  (J  >  0)  the  slope  of  the  curve  in 
Fig.  3(a)  has  the  two  asymptotic  values  indicated  by  the  dashed  lines.  The  ferromagnetic  Curie 
temperature  T^^  is  defined  as  the  position  of  maximum  slope  of  the  x  vs  l/T  curve.  The  anti¬ 
ferromagnetic  N^el  temperature  T^j  is  defined  as  the  position  of  the  maximum  in  the  suscepti¬ 
bility  curve.  With  these  definitions,  the  Curie  and  Neel  temperatures  have  been  determined  as 
a  function  of  x  and  y.  These  results  are  shown  in  Fig.  4. 

The  heat  capacity  is  shown  for  the  same  parameters  in  Fig.  5.  The  positions  of  T^^  and  T^^, 
as  determined  from  the  susceptibility,  are  also  indicated.  The  maximum  of  the  heat  capacity 
curves  differs  somewhat  from  T^^  and  Tj^,  as  determined  from  the  susceptibility. 


TABLE  1  (Continued) 

SPIN,  SYMMETRY  TYPE  AND  ENERGY  EIGENVALUES 
FOR  THE  CUBIC  CLUSTER  OF  EIGHT  SPINS 


Representation 


Eigenvalue  (X  =  E/(—2J)) 


+  2X^(1  +  X  +  y)  +  X(-2  +  7x  +  5y  +  3xy  -  x^) 
+  (—  2  +  2y  +  4x^  +  6xy  —  2x^  =  0 


X  =  —  i  (1  +  3x)  ±  i  [5  —  2x  +  x^  +  4y^  -  8y]'^^ 


X  =-2  +  x-y  ±  (4-  7x-y  +  4x^  +  y^-xyl’^^ 


X=— ■j(3+x)±-j[5— 10x  +  9x^  —  8xy  +  4y^]^'^^ 


X=-2x 


X^  +  X^(5  +  X  +  3y)  +  X(—  i  +  18x  +  14y  —  9x^  +  ixy 
-y^)  3(— 9  +  3x  +  5y  +  3x^  +  6xy  —  y^  —  3x^ 

—  3x^y  +  3xy^  —  y^  =  0 


X=-(l  +  2x)±  [2-2x  +  x^-2y+y^l’/2 


X  =— 3x 

X  =— (1  +  X  +  y) 


X  =  — 2—  X 


Fig.  3.  Susc«ptlbillty  vs  2J/kJ  for  (a)  J  >  0  and  (b)  J  <  0  and  x  =  y  =  0. 
Th*  Curie  and  Ne**l  temperatures  are  indicated.  See  text  for  constant  of 
proportionality,  a. 
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in.  DISCUSSION 


The  calculation  of  the  magrietic:  jaar-opex-ties  for  an  8— spin  cu.t 
esting  features.  The  antlfer*x*om.a|g-ne-tl.c=  state  gives  a  goodl  appr-i 
netic  state  of  the  infinite  solid.  Tlxe;  sxxsosptiloility  and  heat  cap 
high-temperature  limits  as  the  injfin.its  lattice.  For  x  =  y  =  O ,  t 
determined  from  the  susceptihility  of  ttxs  8— spin  cluster  is  given 
with  the  Bethe,  Peierls,  "Weiss  value  of  2. 0-1  (see  Ref.  i).  The  fe 
hand,  is  not  well  represented  hy  the  small  cluster.  The  heat  ca] 
high-temperature  limits,  but  the  susoeptilaility  has  the  correct  1; 
temperature  case  (paramagnetic  regflon.)  .  In  addition,  the  Ourie 
given  by  kT^/j  =  0.90  compared  vritli.  tlaa  ^Bethe,  Reierls,  "Weiss 
In  order  to  investigate  the  nuro."t>ax*  of  spins  required  to  repr« 
of  an  infinite  solid,  the  Heisenberg  approximation  (i.e.,  assume 
carried  out  for  clusters  of  8,  8-4  and  2-18  spins.  The  results  for 

Fig.  6.  The  Curie  tenaperature  for  tlxase  dusters  'was  found  fror 
the  results  are  shown  in  Table  II.  A.  ratHer  slow  convergence  is 
netic  case.  To  put  it  another  -way,  2-18  spins  are  too  fe-w  to  repr 
lattice.  The  results  for  the  cluster  sHo'w  that  fewer  spins  are  ne 
magnetic  state  than  the  ferromagnetic:  state  .  Comparison  -with  tl 
not  useful  for  the  antiferromagnetic  case.  The  Heisenberg  apprc 
antiferromagnetic  state  which  get  worse  as  1-3"  increases;  that  is 


2j/kT  — »- 


Pig.  6.  Heat  capacity  for 
Mei  sen  berg  approximation  ‘ 
21  d  spins. 
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TABLE  II 

VALUES  OF  Tr  CALCULATED 

IN  THE  HEISENBERG  APPROXIMATION 

FOR  CUBIC  CLUSTERS  OF  VARIOUS  NUMBERS  OF  SPINS 

N 

kTc/J 

8 

1.3 

64 

1.8 

216 

2.2 

00 

3.0 

Another  interesting  feature  of  the  8-spin  cubic  cluster  is  the  existence  of  an  antiferromag¬ 
netic  state  for  J  >  0.  For  example,  for  y  =  0  and  x  <  — 1/3,  an  antiferromagnetic  state  exists 

4 

with  J  >  0,  which  corresponds  to  a  "classical"  spiral  state;  that  is,  nearest  neighbors  are 
mostly  parallel  to  take  advantage  of  positive  J  and  next  nearest  neighbors  are  mostly  antipar¬ 
allel.  An  examination  of  the  susceptibility  curve  in  the  region  y  =  0  and  x  ^  — 1/3  shows  that  the 
system  begins  to  order  ferromagnetically  because  of  the  greater  weighting  factor  for  the  S  =  4 

state,  but  at  low  temperatures,  the  system  becomes  antiferromagnetic  because  an  S  =  0,  rT 

7  ^ 

state  is  lower  in  energy.  The  condition  on  x  and  y  for  the  accidental  degeneracy  of  the  S  =  4, 

rj  state  and  an  S  =  0,  state  is  given  by 

3(1  +  4x^  +  3x^)  +  7y^  +  12x^y^  +  2y^  +  3x^^y^  +  Zx^y^  =  0  .  (5) 

In  Fig.  4,  the  dashed  branch  of  the  y  =  0  curve  indicates  this  possibility  of  a  "high"  temperature 
ferromagnetic  state.  The  y  =  ±0.3  calculations  also  allow  this  possibility  but,  for  simplicity, 
these  branches  are  not  plotted  in  Fig.  4. 

It  is  also  possible  to  select  values  of  the  parameters  x  and  y  so  that  an  S  =  0,  state  is 
the  lowest  antiferromagnetic  state  of  the  system  for  J  <  0.  The  classical  analog  to  this  situation 
is  discussed  in  Ref.  4. 
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APPENDIX  A 

EIGENFUNCTIONS  OF  THE  8-SPIN  PROBLEM 


The  notation  is  as  follows: 

spin  configuration  (137)  =  ^(1)  a(2)  /3(3)  a{4)  a(5)  a(6)  /3(7)  a{8) 

in  which  Q!(i)  spin  "i"  up  and  ^(i)  -*  spin  "i"  down.  Only  the  maximum  Mg  state  is  listed,  since 

the  lower  states  can  be  obtained  by  application  of  a  lowering  operator.  For  the  represents- 

tions  r,,  only  one  eigenfunction  is  given;  the  other  can  be  obtained  by  taking  the  complex  con- 
^  ±  ± 

jugate  of  the  given  function.  For  and  Tg,  only  one  function  is  given;  the  other  two  can  be 
obtained  by  performing  the  appropriate  symmetry  operations.  When  the  same  representation 
occurs  several  times,  a  lower  case  Greek  letter  is  used  to  distinguish  the  representation.  An 
"i"  is  used  to  label  the  state  in  the  case  of  a  multidimensional  representation. 

The  states  are  eigenstates  of  and  S  but  not  necessarily  of  X.  The  wave  function  is  a 
sum  of  the  spin  configurations.  The  coefficient  of  each  spin  configuration  is  given  in  Table  A-I. 
The  absolute  value  of  the  square  of  the  wave  function  is  the  last  entry  in  the  table. 
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Spin 

Configuration 

(1234) 

(5678) 

(1256) 

(3478) 

(2367) 

(1458) 

(2347) 

(1348) 

(1245) 

(1236) 

(3678) 

(4578) 

(1568) 

(2567) 

(1368) 

(2457) 

(1278) 

(3456) 

(1467) 

(2358) 

(1357) 

(2468) 
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APPENDIX  B 
MATRIX  ELEMENTS 


In  this  appendix  the  matrix  elements  of  the  form  (r.  |  e  +  xp.  +  yvl  Tj)  are  given.  The  non¬ 
degenerate  states  can  be  obtained  directly  from  Table  I,  so  only  the  matrix  elements  for  degen¬ 
erate  states  are  listed. 

For  S  =  2 


r+a  r+/3 


3  5 

-2  +  2  ^-y 

,, 

—  (1  -x) 

—  (1  -  x) 

-  -1  (1  +  x)  +  y 

rja.i  rjof.ii  rj^.i  Tj/J.ii 


rja,  i 

(1  -x) 

0 

0 

1 

3 

1 

rjof,  ii 

0 

(1-x) 

-  y) 

0 

r’^.i 

0 

-w(l  -  y) 

0 

0 

rj^l.ii 

1 

e 

1 

U 

0 

0 

r+a,i 

r+o,i 

(1  +  x)  -  y 

^(l-y) 

•VT 

Y(i-y) 
_ 1 

-  y  (1  +  x)  +  y 

For  S  =  0 


r+a 


